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c/3 ' The ground-state phase diagram of frustrated S — 1 XX Z spin chains with the competing 

nearest- and next-nearest-neighbor antiferromagnetic couphngs is studied using the infinite- 

S ' system density-matrix renormahzation-group method. We find six different phases, namely, 

', ■ the Haldane, gapped chiral, gapless chiral, double Haldane, Neel, and double Neel (uudd) 

P^ . phases. The gapped and gapless chiral phases are characterized by the spontaneous breaking 

^ , of parity, in which the long-range order parameter is a chirality, Ki=SfSf_^_-^ — Sj Sfj^-^, whereas 

the spin correlation decays either exponentially or algebraically. These chiral ordered phases 
appear in a broad region in the phase diagram for A < 0.95, where A is an exchange- anisotropy 

^ ' parameter. The critical properties of phase transitions are also studied. 

G^ ■ 

in 

^^ 
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§1. Introduction 



'^ ' Frustrated antiferromagnetic quantum spin systems have attracted considerable attention over 

Q , the decades since they exhibit a wealth of fascinating phenomena in their ground states and low- 

lying excitations. In general, frustration supresses antiferromagnetic correlations and the tendency 
i^i ■ towards the Neel order. Classical systems, for example, often show a helical ordered state in 

;h ' their ground states in the presence of strong frustration. In quantum systems, the interplay of 

frustration and quantum fluctuations plays an important role which causes exotic phenomena, e.g., 
a spin-liquid state and a novel type of spontaneous symmetry breaking. 

In this paper, we study a one-dimensional anisotropic spin system with the antiferromagnetic 
nearest-neighbor coupling Ji and the frustrating next-nearest-neighbor coupling J2. The model is 

* E-mail address: hikihara@phys03.phys.sci.kobe-u.ac.jp 



described by the XXZ Hamiltonian 

^ = E ^P E {sfSf^^ + S^Sl^ + A5f 5f+,) I , (1.1) 

p=l K I ) 

where Si is a spin-S* spin operator at site / and A > represents an exchange anisotropy. Hereafter, 
we put j = J2/J1 (i > 0). 

In the classical limit, 5* ^ 00, the system exhibits a magnetic long-range order (LRO) character- 
ized by a wavenumber q. The order parameter is defined by 

m{q) = Wsie'^^, (1.2) 

where L is the total number of spins. While the magnetic LRO is of the Neel-type {q = vr) when j is 
smaller than a critical value, j < 1/4, it becomes of helical-type for j > 1/4 with q = cos~^ (— l/4j). 
It should be noticed that both the time-reversal and parity symmetries are broken in this helical 
ordered state. When the system has an Xy-like anisotropy (A < 1), the helical ordered state 
possesses a two-fold discrete degeneracy according as the helix is either right- or left-handed, in 
addition to a continuous degeneracy associated with the original U{1) symmetry of the XY spin. 
This discrete degeneracy is characterized by mutually opposite signs of the total chirality defined 
by,S) 

0« = tE^'' (1-3) 

Ki = SfSl, - SfSf^, = [Si X Si+,l . 

Note that this chirality is distinct from the scalar chirality of the Heisenberg spin often discussed 
in the literature^ defined by xi = 'S'i-i ■ Si x Sj+i : The chirality O^ changes its sign under the 
parity operation but is invariant under the time-reversal operation while the scalar chirality changes 
its sign under both operations. Since the frustrated classical chain ( |1.1D always has a planar spin 
order, the scalar chirality xi vanishes trivially. 

In the S = 1/2 case, the ground-state phase diagram of the system has been extensively studied 
either numericallyu'DH'® or analytically JSQ' These studies have revealed that the system undergoes 
a phase transition from the spin-liquid phase to the dimer phase at j = jc(A) with increasing j. 
The critical value jc(A) has been estimated to be j'c ~ 0.241 for the Heisenberg case (A = 1)BH'B* 
The spin- liquid phase at j < jc(A) is characterized by gapless excitations and an algebraic decay 
of spin correlations, while the dimer phase at j > jc(A) is characterized by a finite energy gap 
above the doubly degenerate ground states and by the spontaneous breaking of both the parity and 
translational symmetries. The dimer phase is divided into two regions by the Lifshitz line Jl(A) 
according to the nature of the short-range spin correlations, where jl = 0.5 in the Heisenberg 
case.cF Whereas the correlations are of the Neel type and the structure factor S{q) has a maximum 



at q = TT for j < jl(A), the correlations become incommensurate and the maximum of S{q) occurs 
at q = Q < TT for j > jl(A). This incommensurate character is regarded as the vestige of the 
heUcal order in the classical case. It has also been shown both in the XY and Heisenberg cases 
that there is no chiral LRO for arbitrary j, i.e., the chirality O^ vanishes in the thermodynamic 
hmit.0) 

As is well known, on the other hand, the low-energy properties of the S = 1 spin chain without 
frustration are essentially different from those of the S = 1/2 chain. Haldane conjectured that 
integer spin chains are in a phase with a finite energy gap in the excitation spectrum and with an 
exponential decay of the two-spin correlation functions, while half-odd integer spin chains are in 
the spin-liquid phase.EJ'li3) This gapped phase for integer spin chains is now called the Haldane 
phase, and is characterized by the string order parameteiii3 

\ PYn I 



Ostr = 7Eexp Ei^^i r^f- (1-4) 



Haldane's conjecture has been verified either numerically or experimentally, and is now well- 
established. It is therefore natural to expect that the frustrated S = 1 spin chain might also 
exhibit low-energy properties which are substantially different from those of both the S = 1/2 and 
the S = oo classical chains. 

The aim of this paper is to determine the ground-state phase diagram of the model (|1.1| ) with 
S = 1 spins in the j-A plane. The ground-state properties of the model have been studied in some 
limiting cases. In the Heisenberg case (A = 1), Kolezhuk et al. have shown that, at j = 0.7444(6), 
the system undergoes a first-order phase transition from the Haldane phase to the 'double Haldane' 
(DH) phase, in which the next-nearest-neighbor coupling J2 is dominant and the system can be 
regarded as two Haldane subchains coupled with a weak inter-sub chain coupling Ji.llSll3 The 
string order parameter Ostr vanishes discontinuously at the transition. Meanwhile, in the case of 
no frustration j = 0, it has been shown that a phase transition from the Haldane phase to the Neel 
phase occurs at A = Ac -^ 1.18.E3't3't3'tH) The transition in this limit has been verified to belong 
to the universality class of the two-dimensional (2D) Ising model. It is also known that, in the Ising 
limit, A — > CO, there is a phase transition at j = 0.5 between the Neel phase and the 'double Neel' 
(DN) phase, the latter being characterized by the magnetic LRO of 'up-up-down-down' ('uudd') 
type. Systematic studies of the model in the entire j-A plane, however, has been lacking so far. 

In the previous workJl^P we have shown that in the pure XY case (A = 0) there occurs a phase 
transition from the Haldane phase to the 'chiral ordered' phase in which only the chirality exhibits 
a LRO (Ok > 0) with no helical LRO {in{q) = for Vg). Here we note that, in the pure XY and 
Heisenberg cases, the absence of magnetic (spin) LRO has rigorously been proven for any j and for 
general S < oo.E2P We have also found that there are two types of the chiral phase, the 'gapped' 
and 'gapless' chiral phases, in which the spin correlation decays exponentially and algebraically. 



respectively. These chiral phases spontaneously break only the parity symmetry, with preserving 
both the time-reversal and translational symmetries. The phase defined by this types of symmetry 
breaking is a novel one in quantum systems, while it has been studied quite intensively in various 
frustrated classical systems.El[EHlE3'Oy'E3'EZl'E§E3'E3Ei|) It is thus interesting to further investigate 
the properties of the chiral phases in frustrated quantum systems. 

In order to determine the phase diagram, we numerically calculate the correlation functions 
associated with the order parameters, m{q), O,^, and Ogtn for open chains, using the infinite- 
system density-matrix renormalization- group (DMRG) method. Examining the behavior of these 
correlation functions at long distance, we find the following six phases , the Haldane, gapped and 
gapless chiral, DH, Neel, and DN phases, in the region where j > and A > 0. The obtained 
phase diagram has revealed that the gapped and gapless chiral phases exist in a broad region of 
the diagram where < A < 0.95. 

The plan of this paper is as follows. In § |2[ we introduce the spin, chiral, and string correlation 
functions associated with each order parameter. We also explain our numerical method in some 
detail. Results of our numerical calculations are presented in § |3|. The ground-state phase diagram 
in the j-A plane is constructed and the critical properties of the phase transitions between different 
phases are analyzed. Finally, we summarize our results in § ^. 

§2. Correlation Functions and Numerical Method 

In order to determine the phase diagram, we calculate the spin, chiral, and string correlation 
functions defined by, 

C(0 = {SZ^r/2St+r/2) (« = ^, ^), (2-1) 

Cnir) = {K'lo^r/2l^lo+r/2), (2-2) 

/ lo+r/2-1 \ 

Cstr(r) = (5f„.,/2exp i^ ^ 5|pfo+r/2)> (2-3) 

V j=lo~r/2 J 

which are associated with the helical, chiral, and string order parameters, respectively. Here, Iq 
represents the center position of the open chain, i.e., /q = L/2 for even r and Iq = {L + l)/2 for odd 
r, and (• • •) represents the expectation value in the lowest energy state in the subspace of S'^^o^ai ~ ^ 
with even parity. We have checked by exact-diagonalization calculation that the ground state 
of the system indeed belongs to this subspace. We employ the infinite-system DMRG algorithm 
introduced by WhiteJ^I^S The algorithm is accelerated by using the initial Lanczos vector derived 
from the recursion relation proposed by Nishino and Okunishi.E3t3 

Here we note the following two points concerning the properties of the infinite-system DMRG. 
The first is the dependence of the data on the number of kept states, m. Since no analytic form 
to extrapolate the data to m — > cx) has been established so far, in order to obtain precise results 
we have no choice but to increase m until the data converge. Generally speaking, the convergence 



is slow when the energy gap is small, while it is fast when the gap is large. In fact, we have kept 
up to m = 400 states in the vicinity of phase boundaries and in the gapless chiral phase, whereas 
m = 220 turns out to be large enough in the gapped phases. 

Besides, we need to check the L-convergence of the correlations. In the infinite-system DMRG, 
two spins are added in the center of the chain at each step of the algorithm. We continue the 
procedure until the correlation functions calculated around the center become free from the effects 
of open boundaries. In general, more steps are needed if the system has longer correlations. We 
have performed calculations with up to 3500 steps, which corresponds to the system of L = 7000, 
until the L-convergence of the data of the correlation functions was attained. 

§3. Numerical Results 

Examining the r-dependence of the correlation functions given by eqs.(2.1)-(2.3), we find six 
different phases in the phase diagram, namely, the Haldane, gapped chiral, gapless chiral, DH, 
Neel, and DN phases. Figure 1 shows the obtained phase diagram in the j-A plane. The long- 
distance behavior of the correlation functions in each phase is summarized in Table I. In the 
following subsections, we discuss our numerical results at each transition line in more detail. 

3. 1 Transition to the chiral phases 

Let us first consider the phase transition between the Haldane and chiral ordered phases. In our 
previous workJlH) we studied the phase transition in the pure XY case (A = 0), and found that 
there exist two different types of chiral ordered phase. In the gapped chiral phase, the chiral and 
string LRO's coexist and the spin correlations decay exponentially with a finite energy gap, while 
in the gapless chiral phase, the chiral LRO exists and both the spin and string correlations decay 
algebraically with gapless excitations. We also showed that the gapped chiral phase existed in a 
very narrow region of j between the Haldane and the gapless chiral phases. As j increases, two 
transitions, the Haldane-gapped chiral transition at j = id (A = 0) = 0.473 ±0.001 and the gapped 
chiral-gapless chiral transition at j = jc2(A = 0) = 0.490]'^o!oo5' occur successively. 

In the present calculation, we find that the gapped chiral phase continues to exist for A > 0. 
As an example, we show in Fig. 2 the r-dependence of the spin, chiral, and string correlations for 
A = 0.4 and several typical values of j on log-log plots. One can clearly see from the figure, both 
the chiral and string LRO's coexist at j = 0.590, which indicates that the point j = 0.590 belongs 
to the gapped chiral phase. Performing the calculations with increasing j at intervals of 0.002, we 
estimate the transition points for A = 0.4 to be jd = 0.582 ± 0.002 and jc2 = 0.600 ± 0.003. 

Since jd lies so close to ic2) one may wonder if the appearance of the gapped chiral phase might 
be just an artifact of the numerical calculation. We, however, believe that we can refute such 
suspicion. In our infinite-system DMRG, the only source of uncertainty comes from the truncation 
error due to a finite m. As can be seen in Fig. 2, both the chiral correlations in the Haldane phase 



(i^ 0.580) and the string correlations in the gapless chiral phase (j> 0.600) tend to be enhanced 
with increasing m. This means that the estimated transition point jd (jc2) monotonously decreases 
(increases) as m is increased, i.e., the DMRG calculation tends to underestimate the range of j for 
the gapped chiral phase. We can therefore safely conclude that the two transition points locate 
very near, but are distinctly different from each other, and the gapped chiral phase does exist. 
Analyzing the data for several A, we find the gapped chiral phase in the range < A < 0.8. The 
estimated values of jci(A) and jc2(A) are summarized in Table II. At A = 0.9, no region belonging 
to the gapped chiral phase is observed within the calculation with varying j at intervals of 0.002. 
We, however, cannot precisely determine the point beyond which the gapped chiral phase vanishes, 
since we cannot exclude the possibility that the region of the gapped chiral phase is just too narrow 
to detect even at A = 0.9. 

In order to investigate the critical properties of the Haldane-gapped chiral phase transition at 
J = Jci) we analyze the chiral correlation length ^^ in the Haldane phase. The chiral correlation 
length is estimated by fitting the chiral correlation functions in the range 30 < r < 100 to the 
standard form of an exponential decay, 

C,(r) = A exp(-r/C«), (3.1) 

where ^ is a constant. The j-dependence of the estimated chiral correlation length ^^ in the pure 
XY case (A = 0) is shown in Fig. 3 on a log-log plot. It shows a linear behavior near the critical 
point suggesting the standard power-law divergence, ^^ ~ \j — jci\~'^'^- By least-square fitting using 
the three data points near the critical point, we estimate i'^ = 0.9 ± 0.1. The error is rather large 
since the exponent is sensitive to the change of the critical point jd- The estimated value of v^ is 
slightly smaller than, but consistent with, that of the 2D Ising model, i'^ = 1. 

At the gapped chiral-gapless chiral transition at j = jc2, we extract the exponents of algebraic 
decay of the spin correlation function C^{r) in the gapless chiral phase by fitting C^ir) in the range 
10 < r < 100 to the power-law form, 

C^,{r) = BT-^^, (3.2) 

where i? is a constant. Figure 4 shows the j-dependence of the estimated exponent r]x in the pure 
XY case (A = 0). As can be seen in the figure, the exponent rjx increases from r/^. ~ 0.15 at j > 0.8 
to r]x ~ 0.29 at j = 0.490 as j approaches the critical value, j = jc2 = 0-490to;oo5- ^^ estimate rjx 
at the critical point to be rjxij = jc2) = 0.28 ± 0.03. The large error is due to the rather large error 
in the estimated value of jc2- Although the estimated value of rjx at j = jc2 is not inconsistent 
with the value for the Kosterlitz-Thouless transition, r]x = 1/4, more detailed study is needed to 
identify the universality class of the transition between the gapped and gapless chiral phases. 

Next, we consider the transition between the chiral and DH phases at j = jcs- As mentioned 
in § ffl, the system in the DH phase can be regarded as two Haldane subchains with the exchange 



coupling J2 coupled with the weak inter-subchain coupling Ji. Thus it is natural to expect the 
string correlations on each subchain to be the order parameter of this phase. This generalized 
string correlation, however, does not exhibit a LRO in this phase, but decays exponentially .EJ'Ej) 
Consequently, we have no order parameter at hand to characterize the DH phase. We accordingly 
determine the phase boundary between the gapless chiral and DH phases as the line where the 
chiral LRO vanishes. The data of the chiral correlation for j = 1.0 and several typical values of 
A are shown in Fig. 5. It is clearly shown in the figure that the system exhibits a chiral LRO for 
A < 0.68 while C^{r) decays exponentially for A > 0.70. Hence, we estimate the transition point, 
jc3(A = 0.69 lb 0.01) = 1.0. As can be seen in Fig.l, the value of jc3(A) estimated in this way 
becomes larger as A becomes smaller. Here we note that there is a possibility that the gapped 
chiral phase also exists in the vicinity of the line, i.e., the line where the chiral LRO vanishes might 
be different from the line where the energy gap vanishes. It is, however, very difficult to distinguish 
these lines each other, if any, since we have no means to estimate the energy gap accurately enough 
to do this. For that purpose, it may be necessary to introduce another order parameter which 
characterizes the DH phase. 

Since the z-component of the exchange couplings in the Hamiltonian, SfSf,^ reduces the energy 
barrier between the two discretely degenerate states according to the right- and left-handed chirality, 
we can expect naively that the tendency towards chiral ordering is suppressed by the introduction 
of the z-component. This feature can clearly be seen in Fig.l, i.e., the chiral phases appear in the 
narrower range of j as A becomes larger. While the chiral LRO is found for < A < 0.9 in some 
range of j, we observe no chiral LRO for A = 0.95 and 1.0. From this observation, we conclude 
that, between 0.9 < A < 0.95, there is a multi-critical point among the Haldane, gapless chiral, 
and DH phases. 

3.2 Haldane-DH transition 

Here we consider the transition between the Haldane and DH phases at j = jc4- This transition 
has been studied by Kolezhuk et al. in the Heisenberg case (A = l).ot3' From the numerical 
calculation of the string correlation function Cstr(^) by the DMRG method, they estimated the 
transition point, j = jc4 = 0.7444(6), where the string LRO vanished. They also concluded that 
the transition was first-order from the observation of a discontinuous jump of the string order 
parameter from a finite value to zero. Our results are consistent with theirs; we have also observed 
the jump in the string order parameter at the transition, although the estimated transition point 
j = jc4 = 0.747 lb 0.001 is slightly different from theirs. The reason of this small difference may 
be that the number of kept states, m, of our calculation are larger than that of their calculation. 
The first-order nature of the phase transition persists even away from the Heisenberg line (A = 1). 
In Fig. 6, we show the j-dependence of the string correlation at r — > cxd for A = 0.95 and 
1.20. Clearly, the string order parameter disappears discontinuously at the transition point. The 



estimated transition line runs from the multi-critical point to the Neel-DN phase transition point, 
j = 0.5 in the Ising limit A -^ oo, which will be explained in the next subsection. 

3.3 Transition to the magnetic ordered phases 

In this subsection, we discuss the results for the Haldane-Neel transition at j = jcs and the 
region between the DH and DN phases. Some information is available at some special points in 
the phase diagram. In the Ising limit, A — > oo, a phase transition between the Neel and the DN 
phases occurs at j = 0.5. The system exhibits a magnetic LRO in both phases while the nature 
of two-spin correlations changes from the Neel type (j < 0.5) to the 'uudd' type {j > 0.5). In the 
case of j = 0, it has also been shown analyticallyuJ'tZP and numericallj^'cJ^ that there is a phase 
transition between the Haldane and Neel phases at A = Ac ~ 1.18. This transition is continuous 
and belongs to the universality class of the 2D Ising model. In the limit of j — > oo, since the system 
can be regarded as two decoupled chains with the exchange coupling J2, there must be a phase 
transition of the same universality class between the DH and DN phases at A = Ac. 

We estimate numerically the phase boundary between the Haldane and Neel phases by examining 
the existence of the magnetic LRO of the Neel type. As an example, we show in Fig. 7 our numerical 
results of Cg{r) for A = 1.5 and several typical values of j, which suggests that the transition occurs 
at j = jc5(^ = 1-5) ~ 0.1845. As can be seen in Fig. 1, the phase boundary determined in this 
way smoothly connects the transition point in the limiting case, A = Ac for j = 0, to the Neel-DN 
transition point, j = 0.5 for A — > 00. Since the Haldane-Neel transition at j = belongs to the 
universality class of the 2D Ising type, we expect that the transition in the intermediate region also 
exhibits the same critical properties. We analyze the critical behavior of the phase transition in 
the intermediate region at A = 1.5. From the analysis of the order parameter and the correlation 
length associated with the z-component of the spin shown in Figs. 8(a) and 8(b), we obtain the 
critical exponents, /3 = 0.122 it 0.010 and u = 0.96 it 0.08. We also estimate the critical-point decay 
exponent rj = 0.26 it 0.02 from the data of C^{r) shown in Fig. 7. These values of critical exponents 
are in good agreement with the values of the 2D Ising model, f3 = 1/8, rj = 1/4 and u = 1. 

Next we consider the region between the DH and DN phases. As an example, we show in Fig. 
9 the numerical data of Cg{r) for j = 1.0 and several typical values of A on a semi-log plot. In 
the DH phase (A < 1.6), the spin correlation turns out to be incommensurate (note that Cg (r) in 
this phase has been divided by cos{Qr) in the figure) and decays exponentially, while it exhibits 
a magnetic LRO of the 'uudd' type in the DN phase (A > 2.0). In the DH phase, the correlation 
length associated with the incommensurate short-range order becomes longer as A increases. For 
A = 1.7, 1.8, and 1.9, on the other hand, our infinite-system DMRG calculation does not converge, 
i.e., the bond energy at the center of the chain fluctuates irregularly even after thousands of steps 
of the algorithm. This mean that in the region it is impossible to obtain the bulk properties of the 
system within the framework of the infinite-system DMRG method. The origin of this difficulty 



can be attributed to an approximation employed in the algorithm. As noted in § |2|, the size of the 
system treated in the infinite-system DMRG increases by two at every step. Accordingly, the set of 
m-states kept in a step is, strictly speaking, not the optimal one for the block states in the next step. 
In the absence of magnetic LRO or in the presence of commensurate LRO with the periodicity of 
two or four lattice sites (g = vr or vr/2), this error of the algorithm becomes less and less significant 
as the system becomes larger, while in the presence of incommensurate (or highercommensurate) 
magnetic LRO with a wavenumber g = Q ^ vr, 7r/2, the error remains even in the thermodynamic 
limit and causes a crucial problem in the calculation.!^ We therefore consider that the breakdown 
of our calculation in the range 1.65< A<1.95 suggests that the system has an incommensurate 
magnetic LRO in this region. Prom this conjecture, we believe that the point where the magnetic 
LRO of the 'uudd' type starts to appear, j = 1.0 and A ~ 1.95 (see Fig. 9), is a boundary point of 
the DN phase. Connecting the points estimated in this way, we determine the boundary line of the 
DN phase (the dashed line in Fig. 1). In contrast, we cannot precisely determine the boundary line 
of the DH phase nor directly verify the existence of the phase with an incommensurate magnetic 
LRO. 

§4. Summary 

In this paper we have studied the ground-state properties of spin-1 XXZ chains with the an- 
tiferromagnetic nearest-neighbor and the frustrating next-nearest-neighbor couplings. We have 
performed the infinite-system DMRG calculation to numerically compute the spin, chiral, and 
string correlation functions. Analyzing the r-dependence of these correlations, we have determined 
the ground-state phase diagram which turns out to include the six phases, the Haldane, gapped 
and gapless chiral, DH, Neel, and DN phases (Fig.l). The chiral phases appear in a quite broad 
region in the phase diagram for A < 0.95, while the region becomes narrower as A increases. We 
have also confirmed that the gapped chiral phase does exist in a very narrow but finite region be- 
tween the Haldane and gapless chiral phases. The multi-critical point among the Haldane, gapless 
chiral, and DH phases exists between 0.9 < A < 0.95. Furthermore, we have investigated the 
critical properties of the phase transitions. For the Haldane-gapped chiral transition, the chiral 
correlation-length exponent is estimated to be v^ = 0.9 it 0.1, which is not inconsistent with the 
2D Ising value, z^^ = 1- 

The quantum chiral phases which we have discovered in this and previous work are novel ones 
in quantum spin systems. Accordingly, many interesting open problems about these phases still 
remain. As mentioned in § ||, in the S = 1 case only the chirality exhibits a LRO without the helical 
LRO in the chiral phases while both LRO coexist in the helical ordered phase which appears in the 
classical case, S ^ oo. This means that the quantum fluctuation in the S = 1 chain recovers only 
the time-reversal symmetry with the parity symmetry kept broken. In the S = 1/2 case, strong 
quantum fluctuation eventually recovers both symmetries. In order to gain more insight into the 



effect of quantum fluctuation on each symmetry, it is necessary to investigate how the phase diagram 
and the nature of the chiral phases depend on the spin quantum number S, particularly paying 
attention to the possible difference between integer and half-odd integer S. We are now performing 
calculations in the S = 3/2 case. A preliminary calculation suggests that the gapless chiral phase 
appears in this case as well. 
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Table L The long-distance behavior of the spin, chiral, and string correlations in each phase, "expo." and "power" 
means exponential decay and power-law decay, respectively. 

Haldane gapped chiral gapless chiral DH Neel DN 



CI 


expo. 


expo. 


power 


expo. 


LRO 


LRO 


Cf 


expo. 


expo. 


power 


expo. 


expo. 


expo. 


a 


expo. 


LRO 


LRO 


expo. 


expo. 


expo. 


Cstr 


LRO 


LRO 


power 


expo. 


LRO 


LRO 



Table IL The estimated values of jci(A) and jc2(A) for various values of A. 
0.2 0.4 0.6 0.8 



id 0.473 ±0.001 0.527 ±0.002 0.582 ± 0.002 0.642 ± 0.002 0.715 ± 0.003 
jc2 0.490100^5 0.540 ±0.003 0.600 ± 0.003 0.660 ± 0.003 0.720 ± 0.002 



11 



Fig. 1. The ground-state phase diagram. The circles, squares, and diamonds represent the points where the string 
LRO, chiral LRO, and magnetic LRO vanishes, respectively. 



Fig. 2. The r-dependence of the correlation functions for A = 0.4: (a) spin correlation C^(»") divided by the 
oscillating factor cos(Qr); (b) chiral correlation CK(r); (c) string correlation — Cstr(»"). The number of kept states 
is m = 260 (for j = 0.570), m = 300 (for j = 0.580, 0.590, 0.610) and m = 400 (for j = 0.600). To illustrate the 
m-dependence, we also indicate the data with smaller m by crosses for several cases: in figure (b), m — 180, 220 for 
j = 0.570 and m = 220, 260 for j = 0.580; in figure (c), m = 300, 350 for j = 0.600 and m = 260, 300 for j = 0.610. 
In the other cases, the numerical errors of the data are smaller than the symbols. Arrows in the figures represent 
the extrapolated r — oo values. 



Fig. 3. The log-log plot of the j'-dependence of the chiral correlation length for A = {XY case) in the critical 
region of the Haldane phase {j < j'ci). We plot the data setting jci = 0.4729. 



Fig. 4. The j-dependence of the exponent of algebraic decay of CB{r) for A = in the gapless chiral phase 
ij > Jc2 ~ 0.490). The vertical dashed line represents the phase boundary jc2. 
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Fig. 5. The r-dependence of the chiral correlation function for j = 1.0 on a log-log plot. The number of kept states 
is m = 300 (for A = 0.60, 0.65, 0.72), m = 350 (for A = 0.68) and m = 400 (for A = 0.70). For A = 0.68, where 
the 7n-dependence of the data is rather large, we indicate the data with smaller m by crosses {m = 260, 300). For 
the other A, the numerical errors of the data are smaller than the symbols. Arrows in the figures represent the 
extrapolated r — oo values. 



Fig. 6. The strength of the string correlation function at r ^ oo as a function of j for A — 0.95 (solid circles) and 
1.20 (open circles). The vertical dashed lines represent the transition points jc4 for each A. 



Fig. 7. The r-dependence of the spin correlation function C's{r) for A = 1.5 on a log-log plot. The number of kept 
states is m = 220 for all j. The numerical errors of the data are smaller than the symbols. The solid line represents 
the power- law decay of the correlation at j = 0.1845. 



Fig. 8. (a) The strength of the spin correlation function Cs{r) at r — > oo in the Neel phase {j < jcs) as a function 
of j for A — 1.5. We plot the data setting jcs = 0.1845. We extract the exponent l3 by fitting the data to the 
form, Cs{r -^ oo) ~ \j ~ jcsl^ . (b) The correlation length associated with Cf in the Haldane phase {j > jcs) as 
a function of j for A — 1.5. We plot the data setting jcs ~ 0.183. Fitting the data to the form, ^ ^ \j — icsl"^, we 
estimate the exponent v. 



Fig. 9. The r-dependence of the spin correlation function Cs{r) for j — 1.0 on a semi-log plot. For A < 1.6 (the 
DH phase), Cs (r) is divided by the oscillating factor /(r) = cos(Qr), while it is divided by /(r) = exp(i7r[r/2]) 
for A > 2.0 (the DN phase). The number of kept states is m = 180 (for A = 2.0,2.5), and m = 260 (for 
A = 1.2, 1.4, 1.6). The numerical errors of the data are smaller than the symbols. Arrows in the figures represent 
the extrapolated r = oo values. 
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(b) chiral correlation 
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Fig.2(b) T. Hikihara et al. 



(c) string correlation 



10 % 



^ ^^ooooo333iQiDocxmnm: 



1^ 



(((((((((((((««(C 



O 

I 



<w 



10 



-2 



D 



O 



D 




0.570 
0.580 c^ Jci 
0.590 

0.600 ^ Jc2 
0.610 



10 



1 



r 



10' 



Fig. 2(c) T. Hikihara et al. 
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